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We study analytically the Ising model coupled to random lattices in dimension three and 
higher. The family of random lattices we use is generated by the large TV limit of a colored 
tensor model generalizing the two-matrix model for Ising spins on random surfaces. We show 
that, in the continuum limit, the spin system does not exhibit a phase transition at finite 
temperature, in agreement with numerical investigations. Furthermore we outline a general 
. method to study critical behavior in colored tensor models. 
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Cu' I. INTRODUCTION 

Random matrix models provide a statistical theory of random discretized Riemann sur- 
I faces. The amplitudes of the ribbon Feynman graphs of their perturbative expansion support a 1/N 

' expansion Q| (where N is the size of the matrices) indexed by the genus of the surfaces. In the large 

■ N limit the planar graphs corresponding to surfaces of spherical topology dominate. The genus 

i is related (by the Gauss-Bonnet theorem) to the Einstein-Hilbert action on the two-dimensional 

' surfaces and the large parameter appears as the inverse Newton constant. 

00 ' The planar graphs which lead this expansion proliferate exponentially, like K"^, with the number 

, of vertices n. The sum over random lattices of spherical topology is thus convergent for small 

enough coupling constant and the planar free energy is finite. Planar graphs can be counted 
precisely through algebraic equations [1] , as they are related to trees , the universal structures 
behind such equations. When the coupling constant grows, the free energy becomes dominated by 
r> I graphs with a large number of vertices and exhibits a critical behavior. It is in this regime that 

' the system reaches its continuum limit and the critical exponents can be evaluated. 

In the seminal paper followed by 10, 11], Kazakov et al. solved the two-dimensional Ising 



model on random geometries using a two matrix model. Surprisingly this solution turned out to 
be simpler than the one on a fixed lattice. The importance of this work was further enhanced 



by the discovery of the KPZ correspondence 12l4l4i| which relates two-dimensional conformal field 



theories coupled or not to Liouville gravity, and by the introduction of the double scaling limit of 
random matrices which combines all genera isl- 17|. The back and forth navigation between two 



dimensional statistical systems on fixed and random geometries led to early computations of novel 
critical exponents [3]. Their values have been later confirmed by rigorous probabilists, through 
many developments in particular rewarded by the Fields medals of W. Werner and S. Smirnov. The 
transformation of the two dimensional statistical mechanics landscape was so deep that it could be 
considered a change of paradigm. Statistical models on random geometry appear somehow more 
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fundamental and ordinary statistical physics on a fixed lattice as a quenched version in which the 
fluctuating geometry has been frozen. 

It would be highly desirable to generalize such ideas and results to more than two dimensions. 
Random matrices generalize in higher dimensions to random tensors |l9l42ll|. whose perturbative 
expansion performs a sum over random higher dimensional geometries. But until recently the key 
to analytical rather than numerical results, namely the expansion, was missing. 



That situation has changed with the discovery of such a expansion [224241 colored [2 



27l | random tensors. The amplitude of their graphs supports a 1/iV expansion indexed by the degree, 



a positive integer, which plays in higher dimensions the role the genus played in two dimensions. 
The leading order graphs, baptized melonic [2^, triangulate the d-dimensional sphere. They form 
a summable series and map to colored d-ary trees 2^, 2^. When the coupling constant approaches 
its critical value, the free energy exhibits a critical behavior and, like in matrix models, the colored 
tensor models reach their continuum limit dominated by triangulations with an infinite number 
of simplices. The entropy exponent of the melonic series, analogous to the string susceptibility 
Tstring = —1/2 of the 1-matrix model for the pure gravity universal class, is 7meions = 1/2, (28| . 
Colored random tensors are a promising tractable discretization of quantum gravity in three 



and more dimensions and the subject is developing fast [3Ch33I| . The understanding of the leading 
(melonic) order of colored tensor models allows the study of the coupling of statistical systems 
to random geometries in arbitrary dimension. Concretely one needs to develop an algorithm for 
solving tensor models and study the critical behavior of various particular examples. As some of 
the more powerful tools of matrix models, like the reduction to eigenvalues, are absent in tensor 
models, the algorithm we present here relies on combining Schwinger-Dyson (SD) equations with 
a factorization property characteristic of melonic graphs |28|. 

We illustrate our method on the Ising model coupled to random melonic triangulations. It is 
generated by a (colored) two-tensor model, the straightforward generalization of the matrix model 
studied by Kazakov 0] in arbitrary dimensions. It is presented in the Section |TI1 Unlike in two 
dimensions, we find in the Section IIIII that the model does not exhibit a phase transition in the 
continuum limit. The large limit corresponds to the limit of zero bare Newton's constant. As 
such it is only a limiting case of numerical simulations. However, our results are in agreement with 
numerical simulations for small Newton's constant [3j] associated to the (cold) branched-polymer 
phas^. An open question is therefore to understand why and how the progressive freezing of the 
random geometry would let the ordinary Ising phase transition appear. 



II. COLORED RANDOM TENSOR MODELS 
A. Ising Model on Random Surfaces as a Two-Matrix Model 

A single random matrix model in the large limit can be used to describe the behavior of pure 
two dimensional Liouville gravity, with string susceptibility (or entropy) exponent 7 = —1/2 [l|43|]. 
The free energy F expands in the coupling constant g as F Fng"^, where F^ counts the number 

of planar lattices with n vertices (of given valence). The entropy exponent 7 is deduced from the 
asymptotic of F^: if F„ ~ n~^K"', then, for g close to the critical coupling gc, F ~ [gc — g)'^~^ with 
7 = —b + 3. The continuum, large volume, limit dominated by lattices with a infinite number of 
vertices is reached when tuning g to gc- 

In H the critical behavior of the two dimensional Ising model on random planar quadrangula- 
tions is analyzed through a model of two coupled random hermitian matrices X and Y representing 

^ The numerical studies consider triangulations of the sphere which are not necessarily colorable or melonic. This is 
of no concern as many such ensembles are expected to belong to the same universality class. As already shown in 
[28| . the melonic family shares a number of features with the branched polymers. 
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respectively the up and down Ising spins. The free energy of the model is 



F(c, g) = log j dXdY exp Tr - [[/^ + _ 2cXY + + j^Y^] 



(2.1) 



where c 



-2/3J 



, (3 being the inverse temperature, and J > the Ising ferromagnetic coupling. 



The free energy and coupling constant can be solved explicitly in parametrized form as 



F(c,z) 



ff(c, z) 



2g^{z) L2(3z-1)3 



2 Z + 1 



9iz) 

z 



3z-l 



3z 



1 2(3z4 - 3^2 + 1) 



1 z{l-c 
+ -In ^ , , 

2 g{z) 



3 

+ 4 



(1 - 3z) 



c^z + 3c^z^ . 



(2.2) 



The function F{c,g) becomes critical in g when 

dg^ 

dz 







and the random matrix model reaches its continuum limit. The roots of this equation are 



zo 



-i,2 = ^(lTc-V2), 



Z3A = 



^l/2x 



(2.3) 



(2.4) 



In the physical interval of temperature, < c < 1, only the two first roots zq,! ^-re relevant (as the 
others are well separated), and collapse, zo = zi, at c = 1/4. The corresponding critical couplings 
are 



1 

fo(c) = -— + —, gi{c) 



(2.5) 



The first coupling, 50(c), is the critical line (of the large volume limit) at temperature c < 1/4. 
The second coupling, gdc), is the critical line at temperature c > 1/4. At c = 1/4, the two lines 
meet and a phase transition between the low and high temperature phases of the continuous Ising 
model on a random lattice occurs. Evaluating the free energy in the low and high temperature 
phases, F{c,zo{c)) and F(c, zi(c)) one finds that first and second derivatives of F with respect to 
c are continuous at the phase transition, but not the third, hence the transition is third order. 

Note that F„ is the canonical partition function of Ising spins on a random planar quadrangu- 
lation with n vertices. Away from the critical temperature c = 1/4, Fn behaves like n~^K"' with 
b = 7/2, so that the corresponding susceptibility exponent is 7 = —6 + 3 = —1/2, just like pure 
two dimensional Liouville gravity. As the two roots zq and zi meet at the critical temperature. 



the second derivative d g/dz 



1/4 



vanishes and Fn changes its large n behavior to n K"^ with 



b' = 10/3 and corresponding susceptibility exponent 7 = —1/3. 



B. The independent identically distributed 1-tensor model 



ik 



We denote fcj, for i = 0, . . . , D, the D-uple of integers ki = {ka-i, . . . , kio, /cjOj. .^fcjj+i), with 
= 1, . . . , A'^. This N is the size of the tensors and the large limit defined in 23-24| represents 



the limit of infinite size tensors. We set k 



kji. 



Let V- , V'- , with i = 0,...,D, he D + 1 
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couples of complex conjugated tensors with D indices. The independent identically distributed 
(i.i.d.) colored tensor model in dimension D 24, 2^, 27| is defined by the partition function 



e 

S'-'-'-{^, ^) = EEm. + E n + J^nl^. EU^i- (2-6) 

i=0 ^ ^ i=0 ^ i=0 

denotes the sum over all indices ki from 1 to N. 
The partition function of equation (j2.6p is evaluated by colored stranded Feynman graphs 



[25l . l27l |. The colors i of the fields ^*,'(/^* induce important restrictions on the combinatorics of 
stranded graphs. Note that we have two types of vertices, say one of positive (involving ip) and 
one of negative (involving ip) orientation. The lines always join a t/^* to a ip^ and possess a color 
index. The tensor indices kij are preserved along the strands. The amplitude of a graph with 2p 
vertices and T faces is 

A{g) = (AA)PiV-P^^+-^ . (2.7) 



Any Feynman graphs Q of this model is a simplicial pseudo manifold [23] • The colored tensor 
models provide thus a statistical theory of random triangulations in dimensions D, generalizing 
random matrix models. The vertices, edges and faces (closed strands) of the graph represent the 
D, (D — 1) and (D — 2)-simplices. The simplices of the pseudo-manifold are identified by the 
g-bubbles of the graph, for g = 0, . . . , D, i.e. the maximally connected subgraphs made of lines 
with q fixed colors. Bubbles are local objects (in the sense that each of them corresponds to a 
single simplex of the triangulation) and encode the cellular structure of the pseudo manifold. The 
0-bubbles, 1-bubbles and 2-bubbles of a graph are its vertices, lines and faces. 



The amplitude (j2.7p of a graph computes in terms of its degree u}{Q) [2j] as 



A{g) = {X~X)P N^-TD^'^'-^^ , (2.8) 

where p is half the number of vertices of Q. The degree is a positive integer, hence the expansion 
is dominated by graphs of degree oj{G) = 0. Such graphs, dual to triangulations of the D-sphere 
24| in any dimension, are the analogs of the planar graphs of matrix models in higher dimensions. 
They are called melons [i^ and we will restrict our analysis to the melonic sector in the rest of this 
paper. The melonic graphs can be fully characterized at the combinatorial level by a factorization 
property [2^ which we will recall and use in the sequel. 

The free energy Fj;^''^'{XX) has been analyzed in details in [2^. It is analytic for small values 
of the coupling g = XX. When g approaches gc = pji^jyo+Ti the free energy becomes non- analytic 
and the i.i.d. tensor model reaches its continuum, large volume limit. The singular part of the 
free energy is Kjing ~ (s'c — 9) ^ . These results rely on a self-consistency equation for the connected 
melonic 2-point function JJ^-^-'^- and its relation with the free energy 

[/i-i-d- = l + g ([/"•'i) , gOgF'-'''^ = 1 - C/"-^ , (2.9) 



which we will generalize below. 
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C. The random colored two-tensor model 



We turn now to a model with two complex tensors, both colored, {X^,Y^) (and (X*,y*)). We 
consider the partition function 

g-iV^F(xM,c) ^ z{x,x,y,y,c) = j dX dXdY dY e-^(^'^.^^) , 
D 

S{X,X,Y,Y) = Y,Y1 



XlXt+YlYl-cXlYl-cYlXl 



i=0 k 

X ^ _ D 

jv£)(D-i)/4 X] n ^fc, + ]siD{D-i)/4 X] n -^k, 

k k 



y ^ ■ y ^ ■ 

+ jvD{D-i)/4 X n ^4 + ND(D~i)/4^]m, ■ (^-^o) 

k k 

Again denotes the sum over all indices ki from 1 to A^. 

The Feynman graphs of the action p.lOp are identical with the ones of the i.i.d. model, up to 
the fact that a vertex now involves either A or y tensors (both coming with positive and negative 
orientations). Tensor indices kij and colors are preserved along the strands, and the propagator 
between the X and Y sectors is 

Lines in a graph can join a field X (resp. Y) to a field X (resp. Y), with weight 1/(1 — c^), or 
a field X (resp. Y) to Y (resp. X), with weight c/(l — c^). This is the only coupling between 
A-tensors and y-tensors. As in the i.i.d. model, the free energy organizes in powers of 1/A, with 
melonic graphs dominating the large N limit, so that 

F{x,x,y,y,c)= s{g) x^^ x^^yPyfy ^^^^^^ + 0{N-^) , (2.12) 

melons Q 

with (resp. Px, Py and Py) the number of X (resp. X Y and Y) vertices, L is the total number 
of lines, Lxy (resp. Lxy) the number of lines from a A to a y (resp. from a A to a y) and s{Q) 
a symmetry factor. We have L = {D + l){px + Py), Px + Py = Px + Py- 

The A-vertices and y-vertices represent the Ising spins on random lattices. One, say A, mimics 
the random world of the up spins, while the other the random world of the down spins. On a graph 
Q, the total number of vertices n is the number of Ising spins. Through the propagator, two spins 
up or down are coupled with the same weight 1/(1 — c^) and two opposite spins with c/(l — c^), 
from which one can extract the Ising coupling /3J where /3 is the inverse temperature. Taking 
c G [0, 1] corresponds to the ferromagnetic model J > 0. Setting x = x and y = y, the free energy 
of the tensor model equates the grand canonical partition function of the Ising model on random 
lattices, 

Z{fi,f3J,/3h) = ^e""'^ Zn{l3J,l3h) , 

n 

Zn{PJ,Ph) = s{g) J2 ^'^'^^^'''^'^''^^^^''^ ■ (2-13) 

g,2{px+Py)=n Si=±l 
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Zn is the canonical partition function for the Ising system on random lattices with a fixed number 
of spins n, h is the magnetic field and ^ the chemical potential. The relationship between the Ising 
parameters and those of the tensor model is 

III. ANALYSIS OF THE SYSTEM 
A. Method 

The analysis of the critical behavior of our two tensor model is done in three steps 

• We derive SD equations relating the derivatives of the free energy to the connected 2-point 
functions of the tensor model. 

• We evaluate the 2-point functions using the SD equation relating them to the self-energy S 
(IPI amputated 2-point function). In turn, the free energy S writes in the melonic sector 
in terms of the full 2-point functions thanks to the melonic factorization at large iV (the 
defining property of melonic graphs [2^). 

• We send the random lattices to criticality while keeping the temperature c fixed. This 
is the regime where the free energy (the grand canonical partition function of the Ising 
model) is dominated by the graphs with many vertices, hence corresponds to the continuum 
Ising model. We analyze the behavior of the Ising partition function as a function of the 
temperature c. 

B. Derivatives of the grand-canonical free energy 

As the model depends on two complex tensors, there are four types of 2-point functions, noted 
(XX), (yy), (-^y) and (yx) one must deal with. They are necessarily connected and, due to 
the conservation of the indices along the faces of the graph, have the index structure 

{^k^^. = ^n^.n{^X) , (3.15) 

where ^Ui.'pi denotes Il/fc^i '^"i^Pifc ■ ^ correlation functions depend on the variables x,x,y,y,c. 

Proposition 1. The 2-point functions and the derivative of the free energy are related by 

Proof: Using trivial identities of the type 
we get 

= AT- - N-{XX) + cN-{XY) - ^^^(|=T57,( E% n4,) ' 
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n j^i 

that is 

i-(xx> + c(xy> = x — , 

OF 

1- {XX) +c{YX) = X — . (3.19) 

This imphes 

□ 

Also note that 

cdcF = N-^{-c)dclnZ = {-c)((XY) + (YX)) . (3.21) 



If x = x = y = y = the free energy becomes a function of only one variable (denoted by abuse 
of notation F(c,g)) and equates the partition function of an Ising model in zero magnetic field. 
Denoting U = (XX) = (YY) and V = (XY) = (YX) we get 

gdgF{c,g) = 2{1 -U + cV) , cd,F{c,g) = -2cV . (3.22) 



C. Leirge N Melonic Factorization 



The second SD equation we are using is the classical one relating the connected 2-point func- 
tion and the self-energy S. Taking into account that the propagators connecting two self-energy 
insertions S correspond either to X or to Y , we can write it in matrix form as 

_ ({xx)ipi (xy)ip/^ 
lor 2. yijx)^pi (yy)iP7, 

where all 2-point functions carry the same color index. 



(3.23) 



The simplest melonic 2-point graph |28l | is the graph with two vertices connected by D lines. 
Intuitively, more complicated melons are obtained by inserting repeatedly the same motive (i.e. 
two vertices connected by D lines) on lines. In practice this means that the self-energy S in the 
melonic sector factors as the convolution of D connected 2-point functions, one for each color. This 
is the defining property of melons [2^ and is represented in Figure [TJ The key consequence of the 
melonic factorization is that it makes it possible to close the Equation (|3.23p for the full 2-point 
function. Indeed, 

'xx{XX)^ 
^yx{YX)D 

which can be plugged in ()3.23p to yield 

f{XX) {XY)\ r fxx{XX)^ xy{XY) 
[{YX) {YY)) 



xy{XY)^ 
yy{YY)D 



yx{YX)^ yy{YY)^J \ -c^ \c 



1 



1 



(3.24) 



(3.25) 



At zero magnetic field, the matrix equation reduces to two independent equations 



\J - cV - g{U^^^ 



V -cU - gUVi U^^^ + V^~^ 



1 



. 



(3.26) 



Figure 1. Definition of nielonic IPI graplis. Tfie self-energy S factors into connected 2-point functions 
corresponding to the internal lines. 



D. Critical behavior 



The solution for the noninteracting model (g=0) at arbitrary temperature is 



C/(c, 0) = V{c, 0) = , (3.27) 

1 — 1 — 



which is a boundary condition for the system of equations (j3.26p . Multiplying the first equation 
by [/, the second one by V and subtracting we get 



where we chose the determination of the square root verifying the boundary condition. Indeed, 
U = {XX) is positive, while V{c, 0) > for all c > 0, and V vanishes only if either U = 0, or 
c = 0. The square-root in equation (|3.28p . equals V/U and is smaller than one, since 2{U — V) = 
{{X -Y){X - Y)) > 0. Substituting back into (f3:25D we obtain 



4U,,) = (1 - ,U-)^l - j^^^^ - ,U-{l - JJ^^4^)'"' ■ (3.29) 

Although this equation can be solved analytically for D = 3 (as we will show in the next 
section), this cannot be done in arbitrary D, and in general one can only obtain parametrized 
solutions. As the Equation (j3.29p is an algebraic equation it can potentially lead to non trivial 
critical behaviors. We will see below that this is not the case and the partition function and its 
derivatives are analytic functions in the temperature for < c < 1. 

The reader should keep in mind the following scenario. At fixed temperature we slowly turn 
on g. For every c, there exists a critical constant gdc) for which the tensor model reaches its 
thermodynamic limit, the continuous, large volume phase. The free energy of the continuous 
theory is F{c, gdc)). We will see that F{c,gc{c)) is analytic in c for all < c < 1, that is the 
continuous Ising model on a random lattice in dimension larger than three does not undergo a 



phase transition at finite temperature. The numerical simulations of [3j] indicate that the system 
presents just a very long cross over between the totally ordered phase at c = and the totally 
disordered phase c = 1. We will detail the behavior of the Ising model partition function in the 
infinite temperature limit c — >• 1. 
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1. The 3D 
For D = 3 one can cast equation (|3.29p as 



(3.30) 



Introducing the variable z = g^^^U and squaring we obtain a cubic equation for g^^^ 
whose physical solution (the one which becomes g^^^ = z — when c = 0) is 

1/3 



g 



- " ^ - 3z^ ( 6 ) ' ^^H2J=V 2^(1 -.3)3 - (3-32) 



It is difficult to analyze analytically the critical behavior of p.32p for arbitrary z. In the region of 
small z it becomes 

5^/3 (1 _ c2)^ _ (1 + + 2c^)z^ , (3.33) 
which, going back to the U variable, leads to 

1 = (1 - e-W - (1 + + 2e-W' , = - f,^,.Uc'W' ■ 

and critical behaviors for c — )• 1 

hence ~ (1 — c)^/3 and the initial assumption of small z is self-consistent. 

2. In arbitrary dimension 

a. Change of variables. For arbitrary dimensions one can obtain a parametrized solution of 
(j3.29p . To see this we introduce the couple of variables 



As U = {XX) is expected to decrease with c, we get < z < (^j^_^\y/D ■ The physical interval 
of w is readily obtained recalling (j3.28p . thus < w < 1. In the z,w variables we obtain the 
parametrized system of equations 

c = {l-z^)w-z^w^ , g^/^ = z{l- z^){l-w^) . (3.37) 

When w = 0, one has c = 0, and g^^^ = z(l — z^) which is the fundamental equation of the i.i.d. 
model in the melonic sector [i^. When w = 1, g = and z also vanishes so that c = 1. For 
< tt; < 1 we have 

At fixed temperature, the implicit function theorem yields 

/dw\ _ Dz^-^{w + w^) 
VdzJc ~ l-z^ - Dz^w^-^ 



-i?z^-^(t/; + w;^) < , (^) = 1 - 2^ - Dz^w;^-^ > . (3.38) 

> . (3.39) 
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b. Large volume limit and critical curve. For sufficiently small values g the free energy is 
analytic. We now send the lattices to criticality. We obtain a critical curve gdc), on which the 
free energy F{c) = F(c,gc{c)) describes the Ising system in the large volume limit as a function of 
the temperature. 

The critical curve is defined by the vanishing of dg/dU, so that the expansion of g in term of 
U only begins at the second order. In the z parametrization it becomes 



dzJ c 



D 



[z{l - z^){l - w"^)] 
l-z^ - Dz^w^-^ 



z'^^P{w) - z^QH + R{w) 



n 



(3.40) 



with 



P{w) = D{1 + w^){l + w^-^) + (1 - w^){l + D^w 



^2„..D-l\ 



Q{w) = D{1 + w'){l + w''-') + 2(1 



w 



(3.41) 



The roots of the first derivative write 



z^{w) 



D{1 + w^){l + w^"^) + 2(1 - w' 



D 



(1 + 'u;2)2(l + w^-^y - 4(1 - w;2)2^D-i 



2D{1 + w^){l + w^-i) + 2(1 - w;2)(l + DV^^) 



D{1 + ii;2)(l + u;^-i) + 2(1 - u;2) + dJ (1 + w^y{l + w^-^)^ - 4(1 - w^)V-^ 



zi^{w) = zf{w) 



2D{1 + u;2)(i + y^D-i^ + 2(1 - u;2)(l + DV~^) 



(3.42) 



Starting from c = and increasing the temperature, one sees that the physical solution is zilw), 
being the only one such that zf'{0) = jj^ppj- As long as the four roots of the first derivative are 
well separated the Ising spins cannot undergo a second (or higher) order phase transition (the only 
way for the second derivative ^^f to be zero is for two roots of the first derivative to collapse, as 
they do two dimensions where zq = zi at the critical temperature [^). Note that for u; < 1 we 
have 



0<D^ 



(1 - w^){l - w—f + 2^2(1 + w^~^) (1 - w'^){l + w—f + 2u;2(i + ^ 



l + w'^'il + w^'-^f-Ail-w'^yw 



2\2„,,D~l 



(3.43) 



hence the first two roots can never be equal: zi{w) ^ Z2{w). Furthermore, as 



-D 



(1 + w;2)2(l + w^-^Y - 4(1 - u;2)2^D-i < < D(l + w"^) - (1 - w"^) 



< D{1 + w''){l + w^-^) + (1 - w''){l + D'w''-') - 2(1 - w') 



2„..Z)-l^ 



(3.44) 



we have zi{vj) < 1. It follows that a phase transition can occur if and only if zi{w) collapses with 
zsiw) = 0. This happens if and only if P{w)R{w) = and, as P{w) > D, zi{w) = ^ w = 1 and 

0=1. 
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Substituting zi{w) into (j3.37p we obtain the critical curve {c{w), g{w)) parametrized by w 



2DwD{-D - 1 + (D - l)w;2) - 2w{l + D + {D - l)i 



2{l - D^)w''+\1 - w^) 



Dw\l + + w^) (1 - w""-^) + W (1 + w'^D-'^) - 2wD-^ 



(1 + ^2 



9{w) 



^2w(l - + D{1 + w'^){w + w^) - D^JawD+'l (1 - w^f + (1 + lo^)^ {w + w^f^ 



(3.45) 



(2Dw^ ((1 + w2) + D(l - w2)) + 2w (1 + 13 + (-1 + (i _ ^^2) 



^--D 



displayed for some values of D in the Figure m Using = zi{w)/g^/^{w) and = t(;C/('u;), 

one can also plot the parametric curves {U{w),c{w)), iV {w) , c{w)) and {g{w),U{w)). See again 
the Figure [2j 




0.4 0.6 0.8 1.0 




(a)(3(w),c(w;)) 



2.0 2.5 3.0 3.5 4.0 

(h){g{w),U{w)) 





(c)(c/H,cH) (d)(FH,cH) 

Figure 2. Parametric curves on the critical line for 13 = 3, . . . , 8. 



c. The infinite temperature limit. In the neighborhood of w = c = 1, an expansion to second 
order gives (1 — c) ~ "^^^^(1 — w) — -^^(1 — w)"^ . It can be inverted to solve for the parameter w 
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as a function of the temperature c, 

«^ = l-;^(l-^)-^p^(l-^)' + 0((l-c)=^). (3.46) 
One can then expand ah relevant quantities on the critical line around c = 1, 

,D-i n ^^^^ ^-^^ fD + l\ 1 



9cic) - 2--^^^^ (1 - c)-- , UM - V^ic) j ^ , (3.47) 

to lowest order. The derivative of the free energy along the critical curve close to c = 1 is 



dF . , dF dgc OF 

— {c,gc{c)) = —{c,gc{c)) + — —[c,gc{c)) , 

ac oc dc og 

= -2 V,{c) + 2 (1 - U,{c) + cV,{c)) 



(3.48) 



which gives ^ ~ (1 — c)~^. 

Thus the free energy is analytic for all < c < 1. In the infinite temperature limit c — )• 1, 
the physical root zi{w) collapses with z^{w) = and the free energy acquires a non analytic 
behavior. The Ising spin system on random melonic lattices undergoes a phase transition at 
infinite temperature. 
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